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Abstract. Semi-analytical methods, based on Eulerian perturbation theory, are 

a promising tool to follow the time evolution of cosmological perturbations at 
small rcdshifts and at mildly nonlinear scales. All these schemes are based on 
two approximations: the existence of a smoothing scale and the single-stream 
approximation, where velocity dispersion of the dark matter fluid, as well as higher 
moments of the particle distributions, are neglected. Despite being widely recognized, 
these two assumptions are, in principle, incompatible, since any finite smoothing scale 
gives rise to velocity dispersion and higher moments at larger scales. 

Wc describe a new approach to perturbation theory, where the Vlasov and fluid 
equations are derived in presence of a finite coarse-graining scale: this allows a clear 
separation between long and short distance modes and leads to a hybrid approach 
where the former are treated perturbatively and the effect of the latter is encoded 
in external source terms for velocity, velocity dispersion, and all the higher order 
moments, which can be computed from N-body simulations. 

We apply the coarse-grained perturbation theory to the computation of the power 
spectrum and the cross-spectrum between density and velocity dispersion, and compare 
the results with N-body simulations, finding good agreement. 



1. Introduction 

Understanding the statistical properties of matter inhomogeneities in the universe at 
the percent level is one the main goals for cosmology in the near future. A rehable 
description of the evolution of perturbations beyond the linear regime, i.e. at moderately 
high, and high wavenumbers is indeed a very active field of investigation, since future 
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generation of observations, such as high redshift galaxy surveys [H[2l|3llH[5l|6l[7l[8l|9], 
are going to provide information on several aspects of the cosmological model, and 
with their unprecedented accuracy call for a detailed theoretical framework to compare 
with. As an example, the location and amplitude of the Baryon Acoustic Oscillations 
(BAO) in the wavelength range k ~ 0.05 — 0.25 h Mpc~^ are powerful probes of the 
expansion history of the universe and of the properties of the dark energy [TOl. [TT| [T2| [13] . 
Furthermore, the yet unknown absolute neutrino mass scale will be also more efficiently 
constrained (or detected) by comparing data with a theoretically robust determination 
of the power spectrum (PS) in the high k range (for a review, see [E]). 

There are two approaches to deal with nonlinearities. The more established 
one is to use N-body simulations. However, in order to attain the required percent 
accuracy on such large scales very large volumes and high resolutions are needed. The 
resulting limitation in computer time makes it impossible to run simulations over the 
many thousands of cosmologies necessary for grid based or Markov Chain Monte Carlo 
(MCMC) estimation of cosmological parameters. One is therefore forced to develop 
interpolation methods for theoretical predictions, which limits the practical use of this 
approach to "vanilla" type ACDM models and a restricted set of their variants, as it 
was discussed thoroughly in [151 US El [E] . 

The alternative approach is provided by Eulerian perturbation theory (PT), where 
the Euler-Poisson system of equations describing the self-gravitating dark matter fluid 
is solved perturbatively in the matter density fluctuations (for a review, see [19], for 
recent applications, see for instance [20l EH |22l |23l Elj). After its first formulation in 
the nineties, this framework is now experiencing a renewed interest, mainly thanks to the 
work of Crocce and Scoccimarro [2S1 |2S] , who showed that some classes of perturbative 
corrections can be resummed at all orders, enhancing the range of applicability of the 
approach towards lower redshifts and smaller scales. A number of different semi-analytic 
resummation methods have been proposed [271 i2Hl EHl |29l jSOj [3ll [321 [33] and applied to 
the calculation of the matter PS mainly in the BAO range [311 ISSl |35l [361 EZj • Using these 
methods, nonlinear effects have been computed in a variety of non-ACDM cosmologies, 
such as those with massive neutrinos [38], with primordial non-gaussianity of different 
shapes [39|, or with clusterizing dark energy [101 HI]. A quantitative comparison between 
some of these approaches has been presented in [12], where it was shown that in the 
weakly nonlinear regime, most of these approaches attain a percent level accuracy, 
typically better than 2% for ~ 0.1 h Mpc^^ at 2; = 0. 

The crucial question is now how to extend the range of applicability of these 
methods towards larger wavenumbers. The obvious path is to improve the resummation 
schemes in order to include larger classes of perturbative corrections, not yet taken into 
account (for recent attempts along these lines, see for instance [131 [H]). However, even 
if all the perturbative corrections were computed or, equivalently, if the Euler-Poisson 
system was solved exactly, this would not imply that the dynamics of the ensemble of 
self-gravitating dark matter particles would be described with infinite accuracy at all 
scales. Indeed, PT and all the above mentioned resummation schemes suffer from a 
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fundamental limitation, which will strike in at small scales and late times. The point 
is that these approaches describe matter as a fluid, characterized in terms of density 
and velocity perturbations only. This is achieved by truncating the infinite tower of 
coupled differential equations for the moments of the particle distribution, keeping only 
the continuity and Euler equations, complemented with the Poisson equation relating 
the density field and the gravitational potential. All higher moments, in particular 
the velocity dispersion tensor, a*-', are neglected. This picture is fully consistent as 
long as there are no deviations of particle motions from a single coherent flow, the so 
called single stream approximation. Though for non relativistic dark matter particles 
one expects the contribution of such deviations to be sub-leading, yet their inclusion 
leads to corrections to the PS which should be taken into account in view of future 
data sensitivity. This has been recently emphasized in [JH], where the contribution of 
cr'-', which enters the Euler equation, was estimated using numerical simulations and 
quantified to be up to few percent for the velocity divergence PS in the BAO region at 
2; = 0, and slightly smaller for the density field. In a complementary analysis, Valageas 
[I6] estimated that the impact of the single stream approximation on the PS is larger 
than 1 % for A; ~ 0.23 h Mpc^^ at 2; = and for k ~ 0.44 h Mpc~^ at 2; = 1, therefore this 
assumption should be reconsidered if one aims at pushing PT-derived methods beyond 
such scales. 

Another point which is related to the previous discussion is that to obtain a fluid 
description one has to take averages of the spiky particle distribution function over 
volumes typically containing many particles. In N-body simulations, this resolution scale 
can be identified with the grid size over which one interpolates particle positions and 
velocities to obtain the matter density and velocity fields according to some algorithm 
such as, e.g., the Cloud-in-Cell one. In the PT framework, the averaging procedure is also 
instrumental in obtaining a well behaved PT expansion since, the larger the averaging 
volume, the smother the fluid density field. However, in general, the averaging procedure 
introduces a non zero velocity dispersion, even in case at the microscopic level there is no 
crossing of particle trajectories (see also [Sj). Therefore, working in the single stream 
approximation amounts to shrinking the averaging length to zero, and this is the reason 
why this scale does not appear explicitly in PT results. 

The aim of this paper is to reconsider in detail the role of the averaging procedure 
in PT, by keeping the averaging length explicit and defining a perturbative expansion 
for a coarse-grained distribution function and its related moments, including moments 
of the distribution up to second order (velocity dispersion). If one starts from the 
microscopic description and the standard kinetic equation, the smeared quantities are 
found to satisfy again a Vlasov equation but with a source term which contains all the 
information about the behaviour of perturbations at short distances. This source term 
provides the way the smeared quantities are sensitive to fluctuations on scales smaller 
than the chosen cut-off length scale L. We notice in particular that the averaged velocity 
dispersion tensor is sourced by a contribution, due to the product of internal forces and 
velocity fluctuations inside the averaging volume. This is different from what happens 
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for the analogous microscopic o"*-' whose evolution is dictated by a homogenous equation, 
so that if one starts with vanishing initial conditions for it, there is no velocity dispersion 
which can be generated due to dynamics. 

At a practical level, the splitting of perturbations between long-distance modes 
with wavenumber k ~ 27r/L, and small scale contributions, opens the possibility to 
a hybrid approach where the former are treated perturbatively while the effect of the 
latter, which is expected to be to some extent less sensitive to the cosmological model, 
can be measured from a reduced set of N-body simulations. 

The paper is organized as follows. After a brief summary of the usual microscopic 
nonlinear fluid equations in Section |2| we define coarse-grained distribution and the 
corresponding moments in Section |3] and deduce their evolution equations. Section |4] 
contains a similar analysis for the cumulant function, i.e. the moment generating 
functional, showing in a compact way how coarse-graining generates source terms for 
all higher order moments. The Fourier space version of the coarse-grained equations 
are then described in Sections |5| [6} and [7], including a first non-trivial check of the 
consistency of the picture. In Section [8] we develop perturbation theory for the coarse- 
grained moments in presence of the new source terms, which, in Section |9l are computed 



in PT as a first approximation. In Section 10 our results for the coarse-grained PS 
and the cross-correlation between density and velocity dispersion are presented and 
compared with N-body simulations for different values of the coarse-graining scale, L. 



Finally, in Section 11, we give our conclusions and discuss possible lines of developments 
of the approach. 

2. "Microscopic" description 

Our starting point is the microscopic distribution function for the "elementary particles" , 
given by the Klimontovich density in the one-particle phase space (/i space) 

/i^(x, p, r) = ^ (5(x - Xa(r))(5(p - Pa(r)) , (1) 

a 

where the sum runs over all N particles. The particle's coordinate and momenta obey 
the Newtonian equations of motion, 

x„ = — , Pa = -amV0(xQ, r) , (2) 
am 

where dots indicate derivatives w.r.t. conformal time, and the gravitational potential 
satisfies the Poisson equation. 



,2 - . . 47rG 

7 . 



lXa,rj 



(3) 



m ^ 5(xa - X/3) - p , 

with p the background comoving density (which stays constant as the universe expands). 

Taking moments of the Klimontovich distribution function one can define particle 
number density 

ni^(x,r)=y rf^p/K(x,p,r) = ^5(x-Xa(r)), (4) 
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d^P /A-(x,p,r) = ^(5(x-x„(r)) 

am ^ 

From the second order moment, 

am am ^ am am 



am 



(5) 



(6) 



we can extract the velocity dispersion tensor, (t]^(x, r), 

dh f:^-t;i,(x,r) 



= II^(x-x„(r)) 



pj 



<(x,r)j /i^(x,p,r) 
^"^^^-<(x,r)') , (7) 



am ) \ am 

which is non-vanishing only in those points were two or more particle trajectories cross 
(shell-crossing) . 

The conservation of the density ([T]) along particle trajectories in phase space gives 
the dynamical equation 



— + — — - amV> X — 
or ma ax* ap* 



known as the Klimontovich equation. 



/i^(x,p,r) = 



(8) 



3. From particles to fluid 

In any practical use in cosmology we are interested in scales much larger than the 
mean inter-particle distance. Therefore, the relevant distribution function is obtained 
by volume averaging over the resolution scale. 



/(x,p,- 



V 



/i^(x + y,p,r) 



where >V(|z|) is a window function normalized to unity 
1 



V 



d-'yW 



(9) 



(10) 



and rapidly vanishing for > 1. If a volume V ^ typically contains many 
elementary particles, the distribution / is a smooth one, compared to the spiky 
Klimontovich one. The momentum integral of / gives the smoothed number density 
of particles. 



n['x, T 



f/3p/(x,p,r) = -^W 



L 



If the window function >V(|2;|) is the top hat one {i.e. yV{\z\ 



1 for \z\ < 1, and 



W(|-2|) = for \z\ > 1), then the sum in (11) counts the particles contained in a volume 



V centered around x and n(x, r) is the corresponding number density. 
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The first moment of tlie distribution function defines tlie peculiar velocity field, 

1 



w*(x,r) 



1 1 



c^^P— /(x,p,r) 
a m 



Pa 



n(x, t)V ^ am 



E 



X - 



(12) 



By taking again a top hat filter, f*(x, r) is given by the velocity of the center of mass 
(c.o.m.) of the particles contained in V . The second moment 



X, r 



.3 P' 
dp 



am am 



/(x,p,r) 



gives 



1 i^pL^w 



nfx, t) V 



am am 



X — x^ 



L 



(13) 



v'{x, t)v^{x, t) + a-*^'(x, r) , (14) 



where cr*-' (x, r) is the velocity dispersion of the particles in V around the c.o.m. value 



with 



n{x,T) V ^ 



X — Xr 



L 



Pa 



am 



(15) 



(16) 



While the coarse-grained density is just the volume average of the microscopic one. 



n X, r 



1 

V 



d^yW 



ni^(x + y,r) 



(17) 



the coarse-grained velocity is given by the volume average weighted by the microscopic 
density field, and it is therefore different from the filtered density field. 



1 1 
nfx, t) V 



d^yW 



^i^(x + y,r)t;]^(x + y,r) 



^AHx + y,rj 



Moreover, the velocity dispersion obtained from the second moment of /, a*-' is given 
by two contributions. 



a^J(x,r) 



1 



n(x, r) V 



d^yW 



ni^(x + y,r)a]:^(x + y,r) 



1 1 
f2(x, r) V 



d^yW 



r2A'(x + y, t)6v'{x + y, t)6v^{x + y, r) , (19) 



with 



6v\^ + y,r)= v'j,{^ + y, r) - t;'(x, r) , (20) 

where the first contribution is the microscopic one coming from the crossing of particle 
trajectories (see eq. ([T])), while the second one is given by the dispersion of the velocities 
of all the particles contained in the coarse-graining volume, and can be non-zero also in 
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absence of microscopic shell-crossing. If one would iterate the coarse-graining procedure, 
by averaging / over a larger volume, e.g. (2L)^, the new velocity dispersion would be 
given by an expression analogous to eq. (19), where the role of the microscopic velocity 
dispersion would now be played by cx*-^ and the second term would be given by the 
velocity dispersion among the c.o.m. velocities of the volumes contained in the (2L)^ 
one. According to this point of view, the coincidence between non-vanishing velocity 
dispersion and shell crossing is exact only in the L — t- limit, and becomes more and 
more irrelevant as L grows beyond the typical interparticle distance. Indeed, as we will 
see in eq. (27), the macroscopic velocity dispersion is a dynamical quantity, which is 
generated by the product of internal forces and velocity fluctuations inside the averaging 
volume. 

Applying (8]) on ([o]) we derive the equation 

d d _i , 9 
^„,A7^0(x,r) — 

y 



+ 



dr ma dx 
am 



/(x,p,r) 



L 



d 

VL+/0(x + y, r) g-Sfi^ + y, P, r) 



(21) 



where we have split the gravitational potential into a long-wavelength and a short- 
wavelength part, 

0(x + y, r) = 0(x, r) + 50(x + y, r) 

with 

0(x,r) = l|t/=^l/w(|)0( 
and, moreover, we have defined 

"^/(x + y, p, r) = //^(x + y, p, r) - /(x, p, r) 



y,^) 



(22) 



(23) 



(24) 



The term at the r.h.s. of eq. (21) represents the contribution of the short-wavelength 
fluctuations to the evolution of the long- wavelength ones. Neglecting it, one obtains 
the Vlasov equation, which is the starting point for cosmological perturbation theory 
(see for instance [IH]). In other terms, PT and the resummation methods considered 
so far assume the L ^ limit. In this paper, we will work instead with the complete 
equation, and discuss the effect of the r.h.s. on the evolution of the coarse-grained 
quantities derived from /. 

Taking moments of eq. (21) we obtain equations for the coarse-grained density, 
velocity, velocity dispersion, and so on. To keep the notation as compact as possible, in 
the following we will omit the time-dependence, where obvious. The continuity equation, 
as expected, is not modified. 



— r2(x) + -— 



(n(: 



X i; X, 



0. 



(25) 



On the other hand, the equation for the average velocity reads 



dr 



+ 



-V;0(x) -- I d'yW 



ra(x) dx'^ 
n{'x + y" 



n(x)(T^*( 



x)) 



n X 



VL+/0(x + y) 



(26) 
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Compared to the Euler equation, it has two extra terms. The last one at the Lh.s. 
describes the effect of the velocity dispersion around the mean velocity, whereas the 
last one at the r.h.s. accounts for the distribution of matter inside the coarse-graining 
volume V . In other words, while — V^0(x) represents the gravitational force generated 
by the monopole of the mass distribution in V , the second term at the r.h.s. is the 
contribution to the force of the higher mass multipoles. 

Eq. 26 was derived in [3S], where, in order to close the system, different dynamical 



models were considered to compute cx*-^ and the source term at the RHS. On the other 
hand, in ref . |l5] , the cx*-^ tensor was measured from a simulation, whereas the source term 
was not considered. In this paper we take a different approach, by adding the dynamical 
equation for the long-distance quantity ex*-' and by considering different approximations, 
namely PT and N-body simulations, to deal with the short- distance information encoded 
in the source terms. 



Taking the second moment of eq. (21) we get the equation for the macroscopic 
velocity dispersion. 



9r 



d 
y 

L 



d 



d 

dx^ 



v' + - 
n 



1 d 



dx^ 



72 X 



n X 



5v^{^ + y)VUy + '^^'(x + y)Vi J 50(x + y 



(27) 



The equations for the irreducible third moment, w*-' , as those for all the higher order 
ones, can be obtained in a completely analogous way, see also next section for a compact 
treatment of all higher order moments. 

The term at the r.h.s. of the equation above provides a source for the velocity 
dispersion, which is sourceless if one starts from the Vlasov equation instead than 
from eq. (21). It clarifies the "microscopic" origin of velocity dispersions, which, as 
anticipated, are generated by the product of internal forces and velocity fluctuations 
inside the averaging volume. 

Notice that the smeared gravitational potential 0(x, r) is related to the coarse- 
grained density fluctuation, 5(x) = mn(x)/p — 1 by the Poisson equation 

V'0(x,r) = ^•H2^]^5(x), (28) 

while 50(x, r) satisfies the same equation with the density fluctuations corresponding 
to 5r?,(x) at the r.h.s., with 



5?T,(x) 



(i3p5/(x,p) 



As a final remark we take into account the averaged vorticity field, w 
also receives a non-vanishing source term, given by 
1 



?7,(x) 



Vn(x) X Jv = — V(5 X Jv 



(29) 
V X V. It 

(30) 
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where Jv is the second term at the RHS of eq. (26), and it is also turned on by the a - 



dependent term at the LHS. The source term in eq. (30) is suppressed by the gradient of 



the coarse-grained density fluctuation 5 with respect to the source terms for the velocity 
divergence and for the velocity dispersion. Moreover, the effect of the cT^*-dependent 
term was shown to have a sub leading impact on the density and velocity PS's in |3S]. 
Therefore we will neglect vorticity in the rest of this paper, leaving it for a future 
analysis. 



4. Coarse-grained cumulant 



Before proceeding to computations, we discuss the generalization of eqs. (26) and (27) 



to moments of arbitrarily high order. It is useful to consider the Fourier transforms of 
the distribution functions w.r.t. momentum: 

M(x,l,r) = y'rf3pe^^/(x,p,r). (31) 

The microscopic distribution function /(x, p, r) can be thought as the Klimontovich 
one, eq. ([T]), or a smoothed version of it, defined as in eq. ^ but taking the smoothing 
scale L = Ljjy, where Luy is much smaller than any cosmological scale we are interested 
in. In practice, we are assuming that / satisfies the standard Vlasov equation, that is. 



eq. (21) with vanishing RHS. M(x, l,r) is the moment generating functional, since its 
derivatives w.r.t. 1, evaluated at 1 = 0, give 

M(x,l = 0,r) 
. (9M(x, 1, r) 



n x,r 



dU 



n['X,T)v\'X, T 



1=0 



a2M(x,l,r) 



dljdlk 



n(x,r)(t;-'(x, r)t;''(x, r) + o-''''(x, r)) 



1=0 



(32) 

and so on. Taking the logarithm of M(x, 1, r) we get the cumulant generating functional, 
C(x, l,r) = logM(x, l,r), whose derivatives with respect to 1 give the irreducible 
moments, 

= log[n(x,r)] , 



C(x,l = 0,r) 
.aC(x,l,r) 



1=0 



-0 



2 5'C'(x,l,r) 



dljdlk 



(T^'^'(x,r) 



1=0 



(33) 

and so on. 

By single stream regime one refers to the situation where the fluid momentum is 
single valued at each point in space, i.e. particle trajectories do not cross. In this 
regime, the distribution function takes the form 

/(x, k, r) = ^?(x, r) 6nik - P(x, r)) . (34) 
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Inserting this form in (31), we can derive the expression for M, and then for C in the 
single stream regime 

.1-Pfx,r) 



C(x,l,r) 



am 



+ log5'(x,r) 



(35) 



from which we can verify that, being C hnear in 1, all cumulants of order higher than 
one vanish. 

Assuming that the microscopic theory satisfies the Vlasov equation, we can work 
out the corresponding equations for M(x, 1, r) and C(x, 1, r). The equation for the latter 



IS 



dr^^ dV 



dVdx^ 



C(x,l,r) 



.gC(x,l,r) 9C(x,l,r) 



-iFV^0(x,r) , 



(36) 



where we recognize the term responsible for nonlinear couplings in the last term at the 
l.h.s., and the source term at the r.h.s. Since the source term is linear in 1, it only 
appears in the equation for the velocity (the Euler equation). Moreover, the structure 



of the equation - and in particular of the nonlinear term - ensures that the form (35) 



is a fixed point, that is, nonlinear terms in 1, corresponding to velocity-dispersion and 
higher order moments, are not generated if they are not present initially. 

Starting from the coarse-grained distribution function /(x, p, r), see eq. ([9|, we can 
define a coarse-grained moment generating functional. 



M(x, l,r) = / d'^pe'"^- /(x, p,r 



V 



d^yW 



M(x + y,l,r),(37) 



whose derivatives give n, nv^ 



+ a^^), and so on. The relation between the 



cumulants generating functionals is therefore 

1 



C(x,l,r)=log 



V 



d^yW 



exp[C(x + y,l,r)] 



(3J 



From the derivatives of C(x, 1, r) we get the averaged cumulants: 
^(x, 1 = 0, r) = log [n(x,r)] , 



— I 



9C(x,l,r) 



dh 



X, r 



1=0 



,2 5'^(x,l,r) 



dljdlk 



c^^'(x,r). 



1=0 



(39) 



From eq. (35 ) and (38 ) we can determine the expression for the coarse-grained cumulants 



generating functional when the microscopic theory is in the single stream regime, 
C(x,l,r) = 



log 



V 



d?yW 



g{yi + y, r) exp i 



l-P(x + y,r) 



ma 



(40) 
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from which we see that, if the particle momentum P(x, r) is not uniform inside the 
coarse-graining volume V , all orders in 1 are present in (7, and therefore all coarse- 
grained cumulants are non- vanishing, despite the microscopic theory being in the single 
stream regime. 

The evolution equation for C can be derived using that for C, eq.(36), 

. 52 " 



dVdx^ 
1 



-iV W(x,r) + 



V 



C(x,l,r 
d^yW 



.9C(x,l,r) 9C(x,l,r) 



e'^(^+^'^'^Vi+y50(x + y, r; 



where 



5C(x + y, 1, r) ^ C(x + y, 1, r) - C(x, 1, r) 



(41) 



(42) 



One can check that taking derivatives of eq. (41) w.r.t. 1 gives eqs. (25), (26), (27), and 
so on. 



The second term at the r.h.s. of eq. (41) contains all orders in 1 even when the 
microscopic theory is in the single stream regime -provided g{'x, r) and P(x, r) in eq. (35 ) 
are not constant in the coarse-graining volume V - and therefore it represents a source 
term for all cumulants of order greater or equal than one. Notice that, unlike eq. (36), 
the equation for C is not closed, since it requires input from the microscopic dynamics 
encoded in 6C and 6(f) appearing in the source term. The equation for 6C can be 
obtained by subtracting eq. (41) from eq. (36), 



dr dp 



-i V 



5C(x, \,t)-i 



dVdx^ 
a5C(x,l,r) 9C(x,l,r) 

dV dx^ 



— I 



95C(x, 1, r) 95C(x, 1, r) 
dli dxi 
aC'(x,l,r) a5C(x,l,r) 

dx^ 



X, r - 



1 

V 



d^yW 



g<5C(x+y,l,r)y. 



x+y 



x + y,r^ 



(43) 



and, of course, the two equations (41) and (43), together with the Poisson equations 
for and 6(f), form a closed system describing the same physics as eq. (36) (and the 
Poisson equation for (f>). However, the splitting of C in a coarse-grained part, C, and a 
fluctuation 6C suggests to treat them differently, e.g. by using perturbation theory for 
the former and measuring the latter from A^-body simulations. 



5. Fourier space 



In the following, we will deal with the system of equations (25), (26), (27), where we 
will set Uijk = 0. It is convenient to Fourier transform the equations. 

The equations for the coarse-grained density fluctuation {i.e. the Fourier transform 
of ^(x), where uq is the average comoving number density), the velocity divergence 
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(^(x) = V^-y*) and velocity dispersion are 

^(k) + ^(k) + / d'qA - qi - q2) ^ 5(qi)^(q2) = 

^(k) + H^(k) + ^H2n„5(k) - hkja'^{\i) = -Je(k) 

k"^ qi ■ q2 ^ 



12 



(44) 



- d qid q2 ^^(k - qi - q2) 



2 2 ^(qi)^(q2) - M2j^'-'(qi)5(q2) 

^ 1112 



(45) 



a^^ (k) + 2?^ a'^ (k) = - / d^qid''q2 ^^(k - qi - qa) 



^^^7q0a^^(q2) 



+^^7qi) fe^'^(q2) + g{a'^^(q2 



(46) 



where we have neglected vorticity, i.e., we have expressed velocity in terms of its 
divergence, v^{p) = —ip'/p'^ ^(p)- In the equation for the velocity we have also neglected 



terms of order a^^ 5 djS and higher coming from the last term at the l.h.s. of eq. (26) 



This can be motivated since the non-linearity of the coarse-grained quantities 6, 9, a^^ 
can be reduced at will by taking a large enough averaging scale, L. 

The source terms, — Je(k) and — J^-' (k), are given by the Fourier transform of the 



divergence of the second term at the r.h.s. of eq. (26) and of the r.h.s. of eq. (27) 
respectively. 



The symmetric tensor (T'-' (k) can be decomposed as 



a-(k) = I (6- - ^) a(k) + 2 - ^ + hH^) 



(47) 



where 



a{k) = Tr(a(k)) , S(k) = ^a^^(k) , Tr{h{k)) = k^h'^ik) = . (48) 



Neglecting the traceless and transverse tensor h^^ we can derive, from (46), the equations 
for a and S, 



cx(k) + 2U aik) = - 1 d^'qid^'q^ doik - qi - qa) 



X 



1 + 



qi • q2 



(qi ■ q2j 

9 2 



^(qi)^(q2) + 3 



(qi • q2j 



S(k) + 2HS(k) 



d-^qid^q2 5D(k- qi - q2) 



1 



+ 



qi • q2 

2ql 

qi ■ q2 (^ik-i^2f 



(k ■ q2) 

k'^ql 



(k ■ qi)^ _ k ■ qi k ■ q2 qi ■ q2 

k^ql 



,k ■ qi k ■ q2qi ■ q2 (k ■ qi 



2qt 



k'^ql 



k'^qlql 



k'^qf 



^(qi)S(q2) 
(49) 

^(qi)^(q2) 

^(qi)S(q2)| 
(50) 
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where 



J,(k) = Tr{JM) , M^) 



(51) 



6. A first check 

As a first application of the above equations, we can compute the zero-mode {i.e. the 
volume average in physical space) of the trace of the velocity dispersion. Taking the 



average in eq. (49) we get 



9Aao(A) + 2ao(A) = j d'q [-P,^{q) + 2 PeM] 
where we have switched from r to 



Jo(A) 

n 



(52) 



(53) 



A = log a/ , 

and we have defined the volume averages (zero modes) of a and 

6n{k)ao{X) ^ (a(k)) Sn{k)f,{X) ^ (J.(k)) . (54) 

Notice that the brackets indicate a spatial average over the total volume, Vtot = 
{27t)^6dO^), '^ot to be confused with the volume V ^ over which we define the 



coarse-grained quantities. The power spectra in eq. (52) are defined as 



o^(q2)) 



(-^E(q2))=5z,(qi + q2)ns(gi: 



(55) 
(56) 



Consider now the r.h.s. of eq. (52 ) in two opposite limits. If one takes the coarse-graining 



scale L to zero, the source term does not contribute, as it encodes physics inside the 
coarse-graining volume, which shrinks to zero. The remaining contributions are given 
by the cross-correlators between the peculiar velocity field and the velocity dispersions 



a and S, eqs. (55) and (56). However, in the same limit L — )■ 0, the latter receive 



non-vanishing contributions only from the microscopic dispersion velocity defined in 
eq. ([7]), i.e., only from departures from the single stream regime at small scales due to 



the crossing of particle trajectories (see also eq. (19)). This regime is out of reach for 



perturbation theory and semi-analytic resummation methods, and therefore one should 
resort to numerical simulations in order to compute these contributions. 

In the opposite limit, i.e. when L'^ ~ — )■ Vtot 00, the situation is just the 



opposite. The contribution of the integral at the r.h.s. of eq. (52) vanishes, since the 



coarse-grained fields, and their power spectra, vanish for momenta larger than 1/L — )■ 0, 



therefore only the contribution from the source remains. From eqs. (27) and (49) we 
have 



J,(k) = 2 1 W(kL) / d\,d'q2 Soik - qi - qa) 4 Sv^iqi) (50(q2) 



(57) 



where W(kL) is the Fourier transform of the filter function. 
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We can approximate the source (57) by using linear PT, which will turn out to be a 
good approximation in the large L limit (see below). The velocity fluctuation can then 
be expressed in terms of its divergence as 5v^{q) = —iq^/q^ ^^(q)- As for the fluctuation 
of the potential, we will use the Poisson equation which reads 

50(q) = -^^^]^5(q)(l - n;(qL)) . (58) 

Putting all together, and using the property VV(kL) — i- 1 as A; — )■ 0, we obtain the 
zero-mode of the source term as 

q^ 



f,{\) = -Sn'Qm I d'q^^ (1 - W(qL)), (59) 

with 



- 5(q2)) = 5D(qi + c^2)Ps8e{qi) • (60) 



The use of the linear approximation to the PS is justifled by the 1/g^ factor, which 
improves the UV convergence of the integral in (59). In this approximation, and using 
the linear limit of the equations for 5 and (56*, which give 

log 
d\ 



dxPseseiq) = -2(1 + ) Pseseiq) + 3 n^Psseiq) , (61) 



with 

^g(qi) 6e{q2) 

^ n n 



)=6D{qi + q2)Ps9 50{qi), (62) 



eq. (52) can be integrated exactly, to give 

ao(A) J d'q^^^'^il - W(qL)) ^ {Sv\^))nn.tK (for L ^ oo) . 

(63) 

The last equality in the equation above underlines the fact that the result coincides with 
the one obtained by taking the spatial average of the square of the fluctuations of the 
peculiar velocity, as given by linear perturbation theory. In other terms, this solution 
coincides, as it should, with the result which can be derived from the deflnition of the 
macroscopic velocity dispersion in eq. (19), assuming no shell crossing has taken place 
(< = 0). 

We stress that the above computation is a non-trivial consistency check of our 
formalism, which is not passed by the usual way of formulating eulerian perturbation 
theory. Indeed, the latter implicitly assumes a vanishing coarse-graining scale, L — )■ 0. In 
this case, the zero mode of the velocity dispersion cxo is initially zero, and it is sourced 
only by the non-perturbative process of small-scale shell crossing. This is, however, 
inconsistent with a direct computation of which gives the non-zero result of 

eq. (63) already at the linear level. In our approach, by contrast, by changing the 
averaging scale L we can highlight the different contributions to the velocity dispersion, 
namely the microscopic one and the one originated by coarse-graining. 



Coarse- Grained Cosmological Perturbation Theory 



15 



From the above example it is also manifest how the failure of "traditional" eulerian 
PT to reproduce this result is to be ascribed to the fact that it neglects the source terms 
of eqs. (45), (49), and (50), which carry crucial information on the short scale physics. 

Finally, we stress that this dynamical mechanism to generate a non-vanishing 
velocity dispersion can be entirely treated in linear PT, and is therefore completely 
different in nature from the non-perturbative mechanism advocated in ref. [19] to amplify 
an initially small a*-' . 



7. Compact form 



We introduce the four-component field v?a(k, as 

/ 5 \ 



Hf 



a 



(64) 



where 



r/ = log(D+(r)/D+(ri„)), (65) 
and / = d\ogD^/ d\oga, with the linear growth factor. Using the equation for /, 



dlogUf 
dr 



H 



3 

.2T 



/ 



(66) 



eqs. (44), (45), (49), and (50) can be cast in a compact form: 



j d^qid^q25Di^ - Qi - (12)1 abdk, qi, q2)^b{(li,'n)^ci(l2,r]) - KiKv) 



(67) 



where 



/ 1 



^abiv) 



■1 



3 S^TTT, 3 f^JT 

■275- 





V 














3^ - 1 





1 


3^ - 1 



(68) 



The explicit expressions for the non-vanishing components of the vertex functions are 
reported in Appendix A. 

The four-component source term is given by 



/ 











1 


Je 






Hf 






fc2 
W2/-2 






v 


fc2 
^2/2 


■h 


) 



(69) 



The different components of the source /ia(k, 77) in (69) can be read from eqs. (26, 27) 
via eqs. (45, 46), and using the Poisson equation (58) to express the fluctuation of the 
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potential, 50, in terms of the density field. Their explicit expressions can be found in 
Appendix B. 



8. Perturbation theory with sources 



The equation of motion for the coarse-grained variables, eq. (67), can be formally solved 

as 

_. rv r 

V5a(k, rj) = gab{v)vr(M) - / ds Qabiv - s) h(k, s) 



(70) 

where the linear propagator gab{f] — s) is defined as the solution of the equation [26] 

{Sabdr, + ^ab) gbciV - s) = 5ac 5d{V " s) , (71) 

with the retarded initial conditions, gab{,x) — j- 6ab for x — )■ O"*", and gab{x) = for x < 0. 
The explicit expression of the linear propagator can be computed following the lines of 
Appendix A of ref. |32]- Taking the approximation Vim/ P — 1, one finds the analytic 
solution 



9abiv) 



B 



ab 



(72) 



where the three matrices A, B, and C are given by 





/ 


2 




-2 





-4 \ 


1 




-3 




3 





6 


5 




















I 













/ 


/ 










1 \ 

















-1 












1 









V 










1 ) 





/ 3 



B 





V 



-1 \ 

-1 



y 



(73) 



The first two terms in eq. (72) correspond, respectively, to the usual growing and 
decaying modes of the two-component system of density and velocity perturbations. 
These are excited by taking the initial perturbation field, (^^"(k) proportional to the 

vectors 

/ 1 \ 

-3/2 


J 



Ua 



( 1 \ 

1 



V / 



V 



(74) 



respectively. The third term in eq. (72) corresponds to two degenerate modes, both 



decaying as exp(— 2?]), which are excited by initial configurations proportional to any 
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linear combination of the vectors 
/ \ 




(75) 



( ^ ^ 

2 -1 

1 ' ^<^= 

V / V 1 y 

Notice that if the field is initially in the growing mode, <^a"(k) = Ma/(k), then the evolved 
field, (7afc(?7)(^™(k) = ^^/(k), is still in the linear growing mode even if VLm/ P 7^ 1 in 
(pl). 



The formal solution (70) can be expanded perturbatively in powers of the vertex 
function -^ahc, 



^,(k,r/) = 5:v.i")(k,r/) 



n=0 



At zeroth order, we just get the first line of eq. (70) 

rn 



'fa \Kv) = gabivWi^) - / dsgabiv ' s) h{k, s) . 

Jo 



(76) 



(77) 



Inserting this at the r.h.s. of eq. (70) we get 



'^\k,r]) = J (fqid^q2 6D{k-q,i-q2)J^ ds Qabiv - s) e" 'ybcd{k,qi,q2) x 

<^i°nqi,^)^l°Hq2,^), (78) 

and, iterating the procedure, we obtain 

i^'^\k,r]) = J d^qid^q2SD{k - qi - q2) ds Qabiv - s) e" 'jbcd{k,qi,q2) x 
{V^^\qus)^f\q2,s) + ^^^\q,,s)^';l\q2,s)) , (79) 

and so on. 

In the following, we will be interested in computing correlators of the perturbation 
fields, such as the PS 

{cpa{ki,v)M^2, V)) = Sniki + k2)Pab{h; V) , (80) 
or the bispectrum (BS) 

((^a(ki, ?7)<^fe(k2, ?7)<^c(k3, v)) = ^oiki + ks + k3)Babciki, h, h] v) ■ (81) 

By using the perturbative expansions above, the computation of these objects reduces 
to the computation of correlators between zeroth order fields, of the form 

(^i°nki,3i)^r(k2,^2)---^i°Hk^,.J). (82) 

In absence of the sources ha such correlators could be decomposed in products of power 
spectra, under the assumption that the initial field <^a"(k) is gaussian. If the source is 
non- vanishing, since it is generally non-gaussian, this is not possible any more. The only 
guidance to the reduction of the m-point correlators above is, besides the gaussianity 
of (f^^(k), the vanishing of the source and field zero modes (notice that the h^ and /;,4 
components contain a factor with respect to and Js, see (69)) 

(^W(k,r/)) = (^r(k)) = (/^a(k,.)) = 0. (83) 
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In order to compute the PS to 0(7^) we need to take into account the generalized PS, 
BS, and trispectrum (TS), defined respectively as 

(^i°nki,3i)^f)(k2,.2)^f(k3,^3)) = 

^^(ki + k2 + k3)i3°f,^(A;i, ^2, ^3; Si, -52, S3) , 

(<^W(ki, .l)#(k2, S2)^(°)(k3, S3))ri°Hk4, S4))conn. = 

Sniki + k2 + ks + k4)7^°^^(ki, k2, ks, k4; si, S2, S3, S4) , 

(84) 

where the trispectrum, as usual, is defined as the connected part of the four-point 
correlator. In terms of these quantities, the PS (80) is given by 



+ J d^qid^q2 6D{k - qi - ^2) ds gadv - s) 
X I e'-fcde{k, qi, q2)B^ae{k, Qu q2] V, s) + (a ^ b) 
+ 2 r ds'gM{ri-s)e'^'' x 

JO 



+ 



lcef{.k, gi, q2hdgh{k, qi, q2)V^g{qi; s, s')Vff^{q2; s, s') 



ds'gdfi 



s — s le 



Jcdeik, qi, q2)lfgh{qi, ^2, k)V^g{q2] s, s)Vhb{k\ s', r/) + (a ^ b) 

+ J d^qid^q2d^pid^p2 Soik - qi - q2) ^ ds gac{r] - s) 

X ^J^ dsgbdiv - s')e'+''6D{k + pi + P2) x 

'ycef{k,quq2hdgh{k,PuP2)7^fgh{(li,(l2,Pi,P2;s,s,s',s') 

+ I ds'gdfis - s')e'+''(5D(-qi + Pi + P2) x 
.Jo 

7cde{k,qi,q2hfgh{quPi,P2)Tgl^b{Pi,P2,q2,k;s',s',s,r]) + (a 4-> 6) 



^5) 



up to 0(7^) terms. 

Using eq. (77) in (84), we realize that in order to compute the PS in (85) we need 
to evaluate correlators such as 

(k)^r(kO) , {^Ti^)hb{k', s)) , {ha{k, s)hb{k', s')) , (86) 

as well as correlators involving more than three f^*"" fields and h sources. To compute 
the correlators involving only (^*" fields we assume that the initial conditions are taken 
at a redshift high enough that the scales larger than the coarse-graining one, i.e. such 
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that k ^ 1/L, are well in the linear regime, and perturbations are in the growing mode 
identified by the Ua vector in eq. ( |74| ). Therefore, we have 

(^-(k)(^r(k')) = 5D(k + k')P^^{k) = 6o{k + k')W\\cL)P\k)u,u, , (87) 

where P^{k) is the linear density-density PS evaluated at the initial redshift. Notice 
the Vr^(kL) factor, making the only difference between the linear PS's for the averaged 
quantities and those for the not-averaged ones. At 2; = Zin we neglect the quadratic 
contribution 0{6n/n f*) in the integral in (18) with respect to the linear one. As usual, 
if the initial field is non-gaussian, also the initial bispectra, trispectra, and so on should 
be given, otherwise, eq. (87) completely defines its statistics. In this paper we will 
work under this assumption, leaving the discussion of primordial non-gaussianities to 
the future. 

The contribution to the 0(7^) PS from correlators involving only c^*" fields {i.e. 
for ha — )■ 0) can be obtained by setting ^^^(q;^, s') — )■ W'^{qL)P^{q)uaUb and 8^^^ = 
T^hcd = in eq. (85). This gives the usual expression for the 1-loop PS (see for instance 
eq. (162) of ref. p^), in which the linear PS have been multiplied by the UV filtering 
functions W^. Therefore, as expected, the contribution from small scale fluctuations are 
entirely encoded in the sources. In the next section, we will show how to compute the 
sources correlators perturbatively in order to reproduce the full 1-loop result, while, in 
Section 10, we will compute the PS components from a N-body simulation, and compare 
with our results. 



9. Perturbative expansion of sources 



In this section, we will compute the PS for the coarse-grained field, Pat, up to second 
order, i.e. keeping contributions up to 0((P™)^, p«"(5P«"^ (5P*")^), where P^"-{k) = 
iy2(kL)P™(A;) and 5P^"(A;) = (1 - W (kL))'^ P'"" (k) , using the exphcit expression for 
the sources /ia(k, 77) given in Appendix B. 

In order to simplify the computation, we will present explicit formulae for the sharp 
cut-off limit, namely, we will consider the filter function 

W^sharp(A;L) = (1 - kL/{27T)) , (88) 

for which the following property holds, 

H^sharp(kL) (1 - W^sharp(kL)) = . (89) 

However, in Section [10} we will also present results obtained with a smooth cut-off. 
We will compute the different components of PS for the fields by expanding also 



the 6(pa fields appearing in the sources (105), (107) perturbatively. It can be done by 



solving the equations of motion for 6(pa, which can be read out from eqs. (43) and (58) 

{Sabdn + ^ab) Sip\k, 7]) = ha{k, 7]) 



+ e'' / d^qid^q25D{k - qi - q2) 



labc{k, qi, q2) 2 ^b{f\i,v)^^M2, v) 
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(90) 

Notice that the source term ha comes with a positive sign, differently from the case of 



eq. (67) for ipa- In the sharp cut-off limit, the field V9a(k, r^) vanishes for k > 27r/L, 
therefore we are interested in computing the correlators at momenta k < 2tt/L. 



At the desired order, the terms of eq. (85) contributing to the PS are 
Pai>{k;r]) = UaUtP'''{k) 



QaM / dsgM{v-s){vTi^)hci{k',s))/6D{k + k') + ia^b) 
Jo 

+ / (isrfs'5fac(?7 - s)5fM(?7 - s')(/ic(k,s)/id(k', s'))/5D(k + k') 
Jo 

+ y d^qid^q2 dn^k - qi - q2) ds gadv - s) 

X J 2P*"(gi)P*"(g2) r ds'gbdiv - s')e'^''jcef{k,qi,q2)ldgh{k,qi,q2)ueUgUfUh 



+ 



4P^"(A;)P*"(g2) / ds'gdf{s - s')e'+''-^cde{k, qi, q2hf9h{qi, q2, k)ueUgUhUb + {a ^ b) 
Jo 



(91) 



The last three lines give the usual expression for the 1-loop PS in which the initial PS, 
that is P*", is replaced by the one for the coarse-grained fields, P*". Therefore, the 
corresponding momentum integrals are cut-off for g ~ 27r/L. 

The small-scale contribution comes from the correlators involving the sources, 
keeping only terms quadratic in P'^^{k) and 5P'^^{k). The 0{6P^"' ) term is obtained from 



the (hchd) correlator by using the source terms (105)-(107) with 5(pa{q,,ri) = Ua5(p'^^{q, 



It gives exactly the same expression as the fifth line of eq. (91), in which both the 
coarse-grained PS are replaced with 5P*". 

The 0(P'"5P*") contribution comes from the {(p^^ hd) correlator, in which the 5ip 



contained in the first lines of eqs. (105) and (107) have to be computed to the first 



perturbative order according to eq. (90). As a result, the expression for the 1-loop PS 



for the coarse-grained field is given by 
Pn(A;,r7)=P^"(A;) 

+ 2e^"P'''{k) 

+ e^"^ (aP^^^[P*"; k] + W(kLfAP^^^[6P'''; k 

for the 1 — 1 [i.e. density-density) component, where 

k'^TX r°° 



AG[P*"; k] + WikL) ( AG[5P'"; k] + le bal 

6 



(92) 



AG[P; k\ 



252 Jo 



dr P(kr) 



12 



158 + lOOr^ - 42r^ + 



l)=^(2 + 7r2)log 



1 + r 



SP'^'iq) 



(93) 
(94) 
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and 



I9" 



(3r + 7x- lOrx'^f 



[l + r'^ — 2rx)'^ 



(95) 



The 5P™ dependent terms in (92) come from the contributions to eq. (91) containing 
the sources. Since we are working in the sharp cut-off hmit, see eq. (89), the sum of the 
momentum integrals in the two AG-terms and in the two AP*^*^ ones in (92 ) reproduces 
the standard one-loop result for the PS at momenta k < 27t/L. The only difference then 
comes from the 1/6 k'^6a'^ term, which vanishes in the limit of vanishing coarse-graining 
scale, L ^ 0, as (5P™ (and P") in that limit. 

A non-vanishing L gives rise to non-vanishing velocity dispersion. Indeed, the 
cross-correlator between ipi and (ps is given by 

P^s{k,ri) = e'" (p'"(A;)H^(kL)AGi3[5P'"; k] + W^ikL)AP^^^[5P'''; k[ 



(96) 



with 



AG,s[P;k] 



6 



k^Ti 



dr P{kr) 



30 



6r H — -{r 



1)3(1 + r2) log 



1 + r 
1 — r 



(97) 



and 



APi^3^^[P; k]= — / drP{kr) / dxP{kVl + - 2rx) 

7 Jo J-i 

(3r + 7x — 10rx^)(x — r) 



(l + r2-2rx)2 
and similarly for the ipi — (p4^ cross-correlator, 

Puik,ri) = e'" (p'"(fc)W^(kL)AGi4[5P'"; k] + W^ikL)APl'l^[SP'''; k] 



(98) 



(99) 



with 



AGi4[P;A;] 



6 



k^TT 

Ia Jo 



dr P{kr) 



38 + 10r2-6r^ + ^(r2 



1)3(1 + r^)loe 



1 + r 



1 — r 



(100) 



and 



AP*^^[P; k] 



2k^n 



/ drP{kr) / dxP{ky/l + r^ — 2 
Jo J-i 

(3r + 7x — lQrx^){rx — l)x 



2rxy 



(101) 




Figure 1. Snapshots of various quantities obtained from the simulation described in 
the text. The left column is the density field, the central one is the trace of the velocity 
dispersion tensor, and the right one the ratio au/vf. The different lines are obtained 
by taking grids of spacing L = 4, 8, and 16Mpc//i, from top to bottom. 



10. Comparison with simulations 

In order to compare our calculations with the numerical results of N-body simulations 
we perform a cosmological simulation using the public available Tree-Particle Mesh 
GADGET-2 code [50]. The simulation consists of a periodic volume of linear size 512 

comoving Mpc which is evolved from z = 90 to z = 0, this volume is filled with 
512^ dark matter particles and no baryons and cooling or astrophysical processes are 
included. The cosmological reference model is a ACDM with the following parameters: 
f^om = 0.3, Qqa = 0.7, h = 0.7, Us = 0.95 and as = 0.85. The initial conditions are 
generated using the Eisenstein & Hu fitting formula [51]. The particle mesh grid used 
to calculate the long range forces is chosen to be 512^, while the gravitational softening 
is 60 kpc/h in comoving units (this value is constant in redshift). 

At each snapshot we extract the positions and peculiar velocities of the dark matter 
particles and we interpolate using a Cloud-In-Cell algorithm these quantities on a 128^ 
grid. We then compute the mass density, the center of mass velocity and the velocity 
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dispersion tensor for each elementary cube of size (4Mpc//;.) , corresponding to the 



quantities given by eqs. (11), (12) or (18), and (15) or ([19j), respectively. Then we group 
these "elementary" cubes in larger cubes, of size (8Mpc//i)3 or (16Mpc//i)^ and for 
each of these larger cubes we compute the mass density, the c.o.m. velocity and the 
velocity dispersion tensor. 

The decomposition of the velocity dispersion given in eq. ( 19 ) illustrates a difference 
between the elementary cubes of of size (4Mpc//;.)^ and the larger ones. For the 
elementary cubes, the velocity dispersion is entirely given by the contribution of the 
second line, namely by the fact that the different particles inside a cube have different 
velocities, whereas the "microscopic" contribution, at this level, would be given by the 
velocity dispersion of sub-particles inside each particle, which we do not resolve. When 
considering larger cubes, on the other hand, the velocity dispersion associated to each 
elementary cube contributes to the "microscopic" term at the first line of eq. (19), 
whereas the dispersion of the center of mass velocities of the different elementary cubes 
contributes to the second line. The procedure could be iterated ad-libitum, with the total 
velocity dispersion at a given volume size playing the role of the microscopic velocity 
dispersion at larger volumes. 

In Figure [T] we qualitatively show some physical quantities extracted from the z = 
output: each line corresponds to the same slice of the (512 Mpc/h)^ simulation seen at a 
given coarse-graining scale. The first column shows the dark matter density, the second 
the trace of the velocity dispersion tensor and the third shows the ratio a^^/vf. 

Following the first column downwards we verify the obvious expectation that 
increasing the coarse-graining scale the density fluctuation field becomes less nonlinear, 
and therefore PT gets more and more well behaved. 

From the second column we see that velocity dispersion is ubiquitous and, as 
expected, correlated with matter density. However, the relevant quantities to gauge 
if the single stream approximation is well motivated are the ratios /viVj, see for 
instance, eq. (14). From the top panel in the third column we see that, as expected, 
the single stream approximation badly fails inside overdense regions and, at 2; = it is 
marginally acceptable elsewhere (the a^^/vf is everywhere ~ 1). However, increasing the 
coarse-graining scale L, we see that the regions where the single stream approximation 
fails get enlarged, namely, the single stream approximation gets worse at larger scales. 
Therefore, the content of the snapshots in Figure [T] can be qualitatively summarized as 
follows. At small coarse-graining scale the single stream approximation is acceptable 
over a large portion of the volume, but the density field is highly nonlinear and therefore 
PT is unreliable. On the other hand, increasing the coarse-graining scale one gets a more 
linear density field, but velocity dispersion becomes non- negligible at larger scales and 
therefore the PT scheme should be enlarged in order to take it into account. 

In Figure [2] we plot two quantities which can be compared to the analytical results 
presented in the previous section. On the left column we have again the density 
field, already shown in Figure [TJ On the right column we have the "macroscopic" 
contribution to the trace of the velocity dispersion tensor, i.e. the second line of eq. (19). 
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Figure 2. On the left column the same density field as in Figure [T] is plotted. The 
right column represents the macroscopic contribution to velocity dispersion, as defined 
by the second line of eq. (f9). 



Indeed, computing the sources in PT, as we did in Section |9} amounts to neglecting the 
"microscopic" contribution to velocity dispersion (as well as to any other deviation from 
the single stream approximation), therefore the first line of eq. (19) is implicitly taken 
to vanish in this computation. On the other hand, it could be included by measuring 
the sources directly from simulations and then following the general scheme described 
in Section |8| This procedure will be described in detail in a forthcoming paper. 

The right column in Figure |2] shows the correlation between the macroscopic velocity 
dispersion and the density field, which is even more evident than for the total velocity 
dispersion. Moreover, we clearly see that increasing the coarse-graining scale L the 
macroscopic velocity dispersion increases. These results are reproduced by our analytical 
computations, as shown in Figure [3] There, we plot the density-density PS, Pn, and the 
cross-correlator between the density and the trace of the macroscopic velocity dispersion, 
Pi3, for different redshifts, z = 1.5, 0, and coarse-graining scales, L = 8, 16Mpc//i. The 
dots represent the corresponding cross-correlators extracted from the simulations. 

A relevant technical issue is the way one implements the coarse-graining procedure. 
As far as simulations are concerned, the natural way to proceed is, as described above, 
to take larger and larger cubes, that is, to use a cubic cut-off in real space. On the other 



Coarse- Grained Cosmological Perturbation Theory 



25 



hand, when it comes to the analytical computations, it is more practical to work in 
momentum space. Here there is no privileged way to take a cut-off. In our computations, 
we used two different schemes: a sharp cut-off, 

W^sharp(fcL) = e (1 - kL/{2Tx)) , (102) 

and a smooth one 

w (l-T\ o f MkL/2) cosikL/2) \ 

the latter being the Fourier transform of the spherical top-hat filter with radius L/2. 

From Figure |3] we see that, at momentum scales approaching 2tt/L, the difference 
between the two cut-off schemes manifests itself in our results both for Pn and P13. At a 
larger coarse-graining scale, L = 16 Mpc/h, the results obtained with the smooth cut-off 
trace the Pn from simulations better than those obtained with the sharp cut-off, which 
- apart from the |/c^ 6al term at the second line of eq. (92 ) - coincides with the standard 
1-loop computation. The cut-off dependence is, however, strongly alleviated by taking 
ratios of quantities computed in the same cut-off scheme, as we show in Figure |4], where 
the ratio Pu/Pu is plotted. 

Finally, in order to estimate the importance of the dynamics at short scales on the 
coarse-grained quantities, in Figure [s] we plot again Pn as obtained in PT (solid line) 



and from simulations (dots), along with the ^P^^-dependent contributions to eq. (92) 
(dash-dotted line) - i.e. the short- wavelength component of the linear PS - and the 
difference between the simulations and the P*"-dependent contributions to eq. (92) 
(open squares). In other words, comparing the squares and the dash-dotted line one 
can estimate the importance of all the short-distance physics not included in the PT 
result, namely, nonlinear effects beyond 1-loop order and deviations from the single 
stream approximation. 



11. Conclusions 



PT and its improvements by means of resummations are a fundamental tool to compare 
cosmological models to observations and will be even more necessary in the near 
future. Therefore, any investigation of their shortcomings and possible improvements are 
welcome. The single stream approximation and some level of coarse-graining are at the 
basis of these approaches. The impact of these theoretical assumptions on the physical 
results has been investigated only sporadically [52], |l5l HE]. The formulation proposed 
in this paper provides a solid and clear framework to discuss these approximations. 

PT, and its variants, can work only for density fields not too far from linearity; 
at low redshifts, such smooth fields can be achieved only through coarse-graining up 
to a scale in the few Mpc range. This inevitably generates a departure from the single 
stream approximation, which should be taken into account. Our approach allows a neat 
separation of modes: by taking an appropriate coarse-graining scale, the long wavelength 
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Figure 3. Comparison between the perturbative computation described in Section [9] 
and simulations. Pu represents the density field PS, whereas P13 is the cross-correlator 
between density and the macroscopic component of the trace of the velocity dispersion 
tensor. The fields have been rescaled according to eq. (64 1. The continues (dashed) 
lines represent the results obtained with the smooth (sharp) cut-off. The dotted line 
is the standard 1-loop result, while the dots are obtained from N-body simulations. 



modes are in the perturbative realm, and their dynamics feels the effect of short distance 
-non-perturbative- modes via external sources terms. 

In this paper, we have presented a first explicit example of how to deal with these 
sources, namely, to expand them in PT. This provides a check that the "traditional" 
1-loop result of PT is recovered (in the sharp cut-off limit) and shows that -even in a 
fully PT computation - velocity dispersion emerges macroscopically as a consequence of 
coarse-graining, and is therefore not exclusively linked to the microscopic phenomenon 
of shell-crossing. The results obtained in this simple scheme are in (surprising ?) good 
agreement with simulations. 

The next step is to extract and characterize these sources from simulations. It 
will be relevant to see to what extent the short-distance information encoded in these 
sources is cosmology-independent. This would open the path to a fast way to compute 
the nonlinear PS for different cosmological models: treat the cosmology-dependent long 
distance modes in PT, and compute the non-perturbative cosmology-independent part 
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Figure 4. Ratios between P13 and Pn. The contimios (dashed) Unes represent the 
results obtained with the smooth (sharp) cut-off. The dots are obtained from N-body 
simulations. 
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Figure 5. The continuous lines and the dots are the results for Pn with the same 
symbols as in Figure [3j The dash-dotted line and the open squares represent the short- 
distance contribution to the perturbative result and from simulations, respectively. 
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once for all in a accurate N-body simulations. 

An alternative approach would be to pursue further the analogy with the Wilsonian 
idea of coarse-graining, and to formulate a Renormalization Group flow by promoting 
the coarse-graining scale L to the flow parameter: it will be the ideal tool to discuss how 
short-distance features impact on the long-distance physics, and therefore to extract the 
-hopefully few- "relevant" parameters linked to the short- distance modes. 
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13. Appendix A 

We list here the explicit expressions of the vertex functions entering the evolution 
equation for the four component fields <^a(k, r^) defined in Section [T] 



7222(A;,gi,g2) 



Iii2{k,qi,q2) 



k\e - ql - ql) 



72i3(fc,gi,g2) 



Aqlql - -ql- qlf 



7214(^,51,92) 



16g| 

Aqlk^ - Sqfql - Aqj + 3(P - qj - qlf) 



16qi 



7323(^,51,92) 



e{k' + qt + 3ql-2e{qj + 2ql)) 



7423(^,51,52) 



7324(^,51,52) 




32g2g4 



7424(^,51,52) 
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9A;6 + (g2 _ g2)2(3^2 _ ^2) _ ^4(^5^2 ^ ^1^2) + ^2(3^4 + 2g2^| + 3g|) 

32g2g4 

with, finally, 7afec(fc, gi, ^2) = 7ac6(^, ^2, gi)- 
14. Appendix B 

Using eqs. ( [26| ) and ( j27| ) one finds that the three non vanishing components of the 
source vector /^^(k, r^) are the following 

/.2(k, r,) = W{kL) J d'q ^ 5y,i(q, r^) -^,v), (105) 

/i3(k, r^) = 3^6" ly (kL) I ci^^ ^) ^^2(k - q, r^) , 

(106) 

h,{K V) = 3^e^ W^(kL) / d\ V) - q, ^) • 

(107) 
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